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We study the space of scaling fields in the Zn symmetric models with the factorized scattering and propose 
simplest algebraic relations between form factors induced by the action of deformed parafermionic currents. 
The construction gives a new free field representation for form factors of perturbed Virasoro algebra primary 
fields, which are parafermionic algebra descendants. We find exact vacuum expectation values of physically 
important fields and study correlation functions of order and disorder fields in the form factor and CFT 
perturbation approaches. 
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1. Introduction 



The calculation of exact correlation functions is one of 
the most complicated and interesting problems in the 
two dimensional integrable quantum field theory. 

Remarkable solutions of this problem were found for 
the free field models including the Ising model with- 
out a magnetic field and for the conformal field theo- 
ries (CFT) P Ej. One of the difficulties in inves- 
tigating massive integrable field theories with a non- 
trivial interaction is related with the fact that the space 
of scaling fields and, in particular, the space of differ- 
ent solutions of form factor equations, has no clear and 
simple description in terms of dynamical symmetry al- 
gebras like the free fermion algebra in the Ising model 
or Virasoro algebra in CFT. Existence of such infinite- 
dimensional algebras drastically simplifies the problem 
and allows applications of powerful algebraic methods. 
For instance in CFTs the Ward identities generated by 
the holomorphic stress energy tensor T(z) permit to ex- 
press the correlation functions of Virasoro descendant 
fields in terms of the correlators of the primary fields 
and to develop the free field approach for their calcula- 
tion in the minimal models of CFT [2]. 

In massive integrable theories the situation with de- 
scendant fields is more difficult. Perturbing away from 
the conformal invariant point 0] one expects that the 
space of scaling fields is in a one-to-one correspondence 
with that of conformal fields. This statement was sup- 
ported by several results on counting solutions of form 
factor equations El m massive integrable field theo- 
ries which demonstrate appearance of conformal algebra 
characters even off criticality 0. This phenomena may 
indicate existence of some deformed conformal algebra 
action in the space of form factors which would allow to 



extract useful information on descendants from knowl- 
edge on primaries (see also recent works in this direction 

El)- 

In this work we study the structure of the space of 
form factors of scaling fields in the parafermionic CFT 
with the central charge 



c=2(N -l)/(N + 2), N = 2,3, 



(1) 



perturbed by the first thermal operator E3j- I n 
the free field approach [111 IT2"] to form factors we de- 
velop further the idea of the paper |13j on relations 
between different form factors of the theory. Namely, 
we introduce actions of simplest parafermionic modes 
in the space of form factors by exploiting the deformed 
parafermionic symmetry |14| of the underlying lattice 
model We check consistency of our construction 

by numerical studying of correlation functions of the 
theory by using a combination of form factor and con- 
formal perturbation theory approaches ^Hj- To do this 
we also compute new exact VEVs of physically impor- 
tant fields. 

2. Space of states of parafermionic CFT '31 
Conformal field theories with the parafermionic symme- 
try describe self-dual critical points of Zjy ( N = 2, 3, . . .) 
generalizations of the Ising model, the last one 
corresponds to the N — 2 case. In a continuum limit 
order parameters Uk{x) (k = 1, . . . , N — 1) determining 
long range correlations of spins have anomalous dimen- 
sions 

_ k(N-k) 
2dk ~ N(N + 2)- (2) 
Their Zjq charge is equal to k, i.e., under the action of 
the Zn symmetry they transform as 



0~k 



, ,kn „ 
0J <Jk, 



neZ . 



(3) 
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Z/v models allow Kramers- Wannier symmetry and the 
theory also contains the disorder parameters fi k (fc = 
1, . . . , N — 1) with the same dimensions 2d k at the self- 
dual point. Operators fi k transform under the action of 
the dual Zn symmetry as 



(4) 



All other fields naturally separating into families with 
the fixed Zn X Zn charges (k, I) behave as 



kn+ln 



n.n' £ Z . 



The parafermions ip k and ip k generalizing usual 
fermions appear in the OPE of the order and disorder 
fields 

a k (z,z)» k (0,0) = \z\- 4d *z A *[M0) + ■■■}, 
a k (z,z)^+(0,0) = \z\- 4d *z A *[ ^(0) +•••]. (5) 

These currents are holomorphic and generate the in- 
finitely dimensional symmetry due to the conservation 
laws 

d^k = , d,# k = . (6) 

We concentrate on the simplest solution of the associa- 
tivity condition for the operator algebra of currents 



(7) 



which corresponds to the conformal field theory with 
the central charge JIJ and the conformal dimensions 



k{N -k) 
N 



(8) 



The fields ipi = tp (i/j-y = ip) are the basic ones in 
the parafermionic algebra. It will be convenient for 
us to consider as well conjugate currents ijfl = ipN-i 

(V^ = i>N-i)- 

In the conformal models the space of states splits 
naturally into a direct sum of subspaces with the spec- 
ified Zn x Zn charge (k,l) 



{F} = © {F}, 



N > m,m >1-N, 



(9) 



where [m, m] — [k + I, k — 1} , m + m £ 2Z. In these no- 
tations parafermionic currents and order-disorder fields 
belong to the following subspaces 



V> e {^}[ 2 ,0] . ^ G { F l[-2,0] > 



(10) 



-t 



°k e {F} [k k] , n k £ {F} [k _ k] . 



Conformal fields are conveniently classified accord- 
ing to the representations of the parafermionic algebra. 



The action of the parafermionic generators A v (A\) is 
defined by the OPE 



-n T[m.m\ •> 



(ii) 



1p(z)<f>[m,m] = z2 z ™ 

Notice that, if 4>\ m ,m] £ Fw., 
mensions (d, d) then the conformal dimensions of fields 



has the conformal di- 



.4 
At 



v K V\m,m 



] ^ {^}[m+2,m] s 
.in 

£{F} 

[m— 2,m] ) 



(12) 



are {d — v,d). 

Basic fields fi k (a k ) are the primaries of the para- 
fermionic algebra. For instance, the following equations 
hold for n > 



^w+„ = aU 

N N 

Ai-- ■■■ - 3t 



+n+l Mfc - ^i+iL + „ ^fe 



1 Mfe = 1 
0. 



(13) 



All other fields of the model are obtained by the action 
of the currents -0 on [i k {o k ). Thus, the space of states 
of the conformal field theory decomposes into a direct 
sum of irreducible representations of the parafermionic 
algebra 

{F} = ©f-Q 1 \fikU , (14) 

with the disorder operators /i k playing the role of the 
highest weight fields. 

However, from the viewpoint of Virasoro algebra, 
each of the spaces \pk]A expands into a direct sum 
of representations with the highest weight fields 0^ ^. 
These Virasoro primaries are, in general, descendants in 
the parafermionic algebra representation. Indeed, let us 
denote as (ip') fik the field with the minimal conformal 
dimension that can be obtained by /—times application 
of the parafermionic generators to fi k (and respec- 
tively for ip) 



(ip^) 1 [i k — A\,_ 1 _ k A\ l _ 3 _ k ■ ■ ■ A\_ k 



(15) 



Its conformal dimensions are easily computed from J5J. 
Then, up to a normalization, the following relations take 
place 



t) k-2l,-k+2l 



., k , 



(16) 



4> k k+2li _ k _ 2l = (VOW/*. U=o,i,...,N-k. 

For example, the physically important thermal opera- 
tors e k — <t>Q k Q are among these primaries. Recall that 
£fe are local with respect to all fields and have the con- 
formal dimensions 



D k = k(k+l)/(N + 2). 



(17) 
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3. Form factors for perturbed primary fields 
We consider the perturbation of the parafcrmionic con- 
formal field theory by the first energy operator 



-4 — 



CFT 



+ A / d x Ei(x) . 



(18) 



The resulting massive theory is again integrable and Zn 
symmetric. Depending on the sign of A the system is 
in the ordered or disordered phase. We fix A > phase 
where the Zn symmetry is destroyed and vacuum ex- 
pectation values of disorder operators are non-zero. 

Alternatively, this theory can be described also as 
a model with factorized scattering of Zjy charged par- 
ticles |SJ. The particles a G {1, . . . ,N — 1} in the Zm 
(N = 2.3,4,.. .) symmetric models have masses 



M a = M 



sm(Tra/N) 
sm(n/N) 



(19) 



Note, that the antiparticle is identified with the parti- 
cle N — a. The scattering matrix of the lightest particles 
a = 1 has the form 



SnW) 



sinh(f + f ) 
sinh(| 



(20) 



N . 



The S matrices for higher particles are also diagonal 
and can be extracted from Sn according to a standard 
bootstrap prescription HJ. 

We follow the algebraic approach [TT1 El H3J HZ] to 
the form factors originated from the analysis of lattice 
ABF model |15|. This lattice model falls into the same 
universality class as the Zm Ising model. In the cor- 
ner transfer matrix approach its dynamical symmetry 
is a deformation of a parafermionic symmetry. The de- 
formed currents 'J (a), as well as deformed primaries 
^m( a ) (which are primary for deformed Virasoro al- 
gebra [12 ) act in the corner transfer matrix space of 
states. Note, that there is an essential difference be- 
tween the symmetries in off-critical situation and in the 
CFT. Thus, there exists only one chiral algebra, the pa- 
rameter a is related with the spectral parameter but not 
with the distance, etc. However, from the mathemat- 
ical viewpoint, the analysis of ABF models, including 
fusions rules, Zn symmetry, the structure of irreducible 
representations etc., is very similar to the conformal sit- 
uation ©-(^UJ). Respectively, as in CFT, the problem of 
the calculation of physical quantities reduces to study- 
ing of operators acting in the direct sum of representa- 
tions of the deformed parafermionic algebra. 

In particular, the form factors are matrix elements of 
local operators in the basis of asymptotic states formed 
by operators diagonalizing the Hamiltonian. A pecu- 
liar feature of the Zjq models is that the Hamiltonian is 



diagonalized by the symmetry generators '5 (^S^) them- 
selves. At the scaling limit the currents become particle 
creation operators B {B^) which in the free fields ap- 
proach 1 1 1 1 can be described as following. 



B(/3) 



BH(3) 



== V ae^ +p ^ a Z a (f3), (21) 

l n -li- i— -* 



V 2s[li 7? „t±i 



rf'2 



ae^- p i+ p *-V a Zl((3). 



These linear combinations contain the exponentials of 
free bosons Z±(f3) and Z±(f3) and the "zero modes" 
operators Pi,i- The operators Z a {(3) and Z\{(5) are as- 
sumed to satisfy the Wick theorem 

({Z ai (A) • • • Z an {(3 n )Zl (#)... Z\ m (O) (22) 
= [] iZ ai (A)Z aj (&)» J] 1 Z l (#K (#)» >< 



%<3 



i<3 



n<<^;(&)-4 (#)>>, 
with the following contraction rules (where /3 = /3i — /?2 ) : 

{{zmzm)) = {{zUMzUfo))) (23) 

sinh(f + &(o - b)) 



sinh ^ 



((ZMztiPz))) = {{zUMZ-bim 

= C t (/3)cosh(^-^(a + 6)). 

The functions CHP) an d their basic properties are 
given in the appendix A. 

The perturbed primary fields 4>m -m corresponds to 
the projectors to the states (($^. TO | and 1$^))- We treat 
the last ones as states created by the deformed operators 
&± m from the vacuum vector 



K»= lim <(a)|0», 

ot— >oo 

«<I>*J= lim «0|*i ro (a). 



(24) 



Let us note that these bra and ket vectors in the corner 
transfer matrix picture correspond to the Virasoro alge- 
bra states created by left and right chiral parts of fields 
-m- The action of the zero modes on the vectors 
(|23} reads as 



P 1 \<Z> k n ))=k\$ k J), P 2 \$ k J) = -m\$ k J) , (25) 

and we assume that the operators Pj commute with 
operators Z a (f3), Z\((3). The normalization of vectors 
1$^)) is chosen according to the prescription 



(<*-J*£) 



m. — ral ? 



(26) 
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which would correspond to the conformal normalization 
of scaling fields 



W*Mo)> = \x 



x^O. 



(27) 



Then the form factors of the Virasoro algebra primary 
fields _ m , which are local with respect to the fields 
jUfc, can be compactly written as following 

«-mK0}> {/3'}>(n,n) = 

n n 

= (^ k . m \l[B(0 j )l[B\(3' j )\^ m }). (28) 

We implied here a shorthand notation for the state with 
n' lightest particles 1 and n antiparticles it 

{/3'})(n,n') = \@U ■ ■ ■ , Pn, Pi, ■ ■ ■ , #t')it...it,i...i 

The symbols of ordered products of particle creation 
operators stand for 

n 

\{B{P j )=B{(3 1 )---B{p n ). 

The equation 128(1 is a natural generalization of the ther- 
mal operators form factors to the m ^ case. With 
the definitions l)21 [l -l)22 [l it has the conventional form 



(<t> k m ,-m\{P},{P'})(n,n) = 



(-If 



■E(Pi-P'j)- 



{ijfij} j 



(2 sinf )™ 

gW((' C + 1 )( a J- b 3)^ m ( a 3+ b j)) X 



((^(^•••^jO^-ruK))- (29) 

The remarkable property of our theory is that form 
factors of perturbed Virasoro primaries ^ includ- 
ing H28fl can be obtained form factors of fields fif. by 
the equation analogous to l(16|l . Indeed, for large a 
the operators $± m (a) from (|24|l behave in the confor- 
mal limit, roughly speaking, as chiral parts of primaries 

-m( e_a > e ~ a )' Respectively, the operators \&(a) in 
CFT limit become the parafcrmionic currents tp(e~ a ). 
Basing on this correspondence, we introduce the analogs 
of rules (|ll|) - (|l(jfl for the action of the simplest Fourier 
modes of deformed currents iff (\ft) as 

lt 2 )) = I^L<»= I™ e-^^t^)!^)), 
«** +2 l = ((Ai±™<\ = lim e^«((^|vl/(a), 

n a — ►— oo 

((C 2 l = PU<h lim e ^((^|*t(a). 

jv a — oo 

Since the operators £> (fit) are, essentially, the scaling 
limits of the deformed parafermions \P (^t), we put 



(Al#><tk,n\{j3},{P}kn,n>)= (30) 
n n' 

lim e -^((^ine(A-)n st ^) 

a^oo A A - 1 - J 

{/3}) ( „, n0 = (31) 

n n 

^lim^ e^ a («|B(a) [] II st (#)K» ■ 



Starting from 128|l we have the i?/v neutrality condition 
2(n' — n) = to + to + 2. For Fourier modes , we obtain 

(4- m ^|{/3},{/nW')= ( 32 ) 

iV 

n n 

Km e -^"<«inm)n 6t (^) 6t wi<)), 

(4UiJM,{/3}> M) = (33) 

iV 

n n 

lim e 1 ^" «^|0t («) [] [] gt ( ^) |<j>* 

CK — > — OO A A A A 



-mil. 



with the selection rule 2(n' — n) = m + m — 2. 

These actions provide maps between solutions of 
form factor equations for fields with the different Zjq 
charges. Note, that the equation (|16fl holds for corre- 
sponding multi-particles form factors, as it can be di- 
rectly verified from l|28|) and (|30I33|) . In this way we 
extend free field formulae for primaries 4>m ^ . 

4- Short distance expansions and VEVs 
As an application for new form factor formulae 1|28[) . 
(|3UI33f) and for checking the consistency of the proposed 
parafermionic actions we analyze the correlation func- 
tions 

G+(x) = (a 1 (x)a+(0)) , G_(x) = ( Ml (x) M +(0)) , (34) 

computing it in the form factor approach (long distance 
expansion) from one side and in the conformal pertur- 
bation theory method (short distance expansion) 

G ± (x)^r- id ^Ci'(r)(M0)), r=\x\, (35) 

from the other side . The structure functions C± (r) 
above allow expansions into perturbation series 

C£ (r) = (C£ (0) + Xr^-^C^ +•••), 

where the coefficients C-j-'m can expressed through 
the integrals of correlation functions in CFT. 

In the first order approximation the main contribu- 
tions to (|35|l come from the primaries J, £i,£2 and the 
W-algebra descendant field E x = W^W^E! 
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ci = i± 

C± = =F 



A7rr 2(l-2n) ( 7 ( m ) 7 (3 m ))1/2 7 ( 4m ) 

2(1-4u)V(2m) 
r 4tl ( 7 (u) 7 (3u)) 1 / 2 Avrr 2 M 2 7(4u) 



(36) 



2 7 (2w) 



4 7 4 (2u) 7 4 (l - u) 



~is2 = -j- 12u A7rr 2 ( 1 - 2 ") M VH( 7 (5 M )) 1 / 2 
y± + 12(l + 2u) 2 7 2 (2ii)( 7 (3w)) 1 /2 



C^ 1 = - r 2 + 4u 



u(l-4u)(7(u)7(3u)) 1 / 2 
(1 + 2u){\ - 2u)j(2u) 

Here and later we denote 

T(z) 1 



700 



T{l-z) 



N + 2 



(37) 



Let us recall that the mass M of the lightest particle 
and the coupling constant A can be expressed via each 
other. The exact relation between them in the confor- 
mal normalization (ei(x)ei(O)) = |.t|~ 4£>1 is known [T%| 
to have a form 



(2^A) 2 = k 4(1 - 2u) 7 (u) 7 (3u) , 

r r(2/jv)r(i - i/n) 



K = M- 



T(l/N) 



(38) 



(39) 



Finally, to compare the CFT perturbation result H35JI - 
1361) with the form factor predictions we have to com- 
pute vacuum expectation values of operators |19| in the 
conformal normalization (|27fl . For \x k these important 
quantities can be derived either from the reflection rela- 
tions |201 |2"T] or via the deformed vertex operators ^7] 



(fl k ) s = U Sk K 2d « X 

dt ( sinh kut smh{N — k)ut 



(40) 



exp / — 



t \ sinh t tanh Nut 



- 2d k e- 2t ) ) . 



Here the index s = 0, N — 1 enumerates the vacuum 
states in the Zjv broken phase. (Note, that the correla- 
tion functions G± do not depend on it.) 

For the general case, _ m — (t/j^) l (ip) '/ifc, where 
m = k — 21, the vacuum expectation values of fields are 
elegantly expressed in terms of (^fc) s as 



C-Jb ( K \ ^ ^_ 2s (t + 1) 7 (^) 



\NuJ 



it 



(41) 



In particular, VEVs of thermal fields St — (tfj) k (ip) k u, 2 k 
correspond to the m — case in the equation above: 



(£fc) s = 



k 2D » (fc!) 2 fj((2k + l)u) 



(Nu) k (2k)\ 



7 2 (^) 



7(m) 



(42) 



n 



\l(§M(2i + l)u) 



Using the method developed in [23 we a ls° obtained 
the exact result for the expectation values of the nor- 
malized descendent fields E k 



(E, 



k s 



(iv + 2) 2 r 2 (i 



N 



r 2 (i 



N ' 



27V r 2 (i-*±i)r 2 (i + ^) 



r(i 



2fc -)r(i ■ 2 



N ■ 



N 



r(2 



M)r(2- 



2fc+2 ' 

N > 



(43) 



Collecting these results we can analyze the short 
and long distance behaviors of correlators. Figs. 1- 
2 demonstrate typical results obtained by numerical 
computations. The lower lines correspond to the short 
distance asymptotic from the CFT perturbation the- 
ory (|35[l - il36|) and the upper lines give the form factor 
decomposition up to two particle. We see that both 
expansions are in a perfect agreement at the interme- 
diate distances and therefore provide good numerical 
data for G± at all scales. This confirms the proposed 
identifications of form factors and also gives one of 
few examples of unitary theories where the correlators 
can be effectively studied in the scheme of Ref. 16 . 

2 . 15 r 
2.1- 
2 . 05 - 



0.5 



1 . 5 



1.95 ; 

1.9- 
1.85- 
1.8- 



Fig 1. Disorder correlator G_(r). N=7 Model. 




0.25 0.5 0.75 1 1.25 1.5 

Fig 2. Order correlator G+(r). N=7 Model. 

5. Discussions 
Algebraic relations among form factors of Zn mod- 
els can be extended for higher Fourier modes of 
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parafermions as well. Though, in general, it is not clear 
how to confirm the proposals, we were able to provide an 
additional check supporting the validity of our construc- 
tion. Namely, we verified the form factor counterpart of 
the quantum equations of motion, which look like: 



C09) 



= (1 - £>i)A ty)ei 
diP = {l-D 1 )\ 



(44) 



The form factors of operators (tp)ei, (4>)e\ are defined 
as in Eqs. (|30[1 - H33[1 . while to compute form factors of 
currents in the left hand side, we expanded correspond- 
ing matrix elements from the lattice theory |14| up to 
first non vanishing terms in the parameter of the lattice 
spacing. Equations l|44|) explain the relations between 
form factors of ip and £\ originally found and discussed 
by F. Smirnov |13j for the Z% model. 

Notice, that for small number of particles the form 
factors obtained from Eqs. (2EJ, (j20Il-l|S3), ffij agree 
up to normalization factors with the results for Z^ case 
[T3) (see also Ref. ££5]) and for case [21] . 
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A. Properties of two particle form factors 
Let the contour C goes from infinity above the real axe, 
then around zero and then to infinity below the real axe. 
Introduce the notation 



Si{x) = exp- 



1 f dt sinh(a; — 27r)t 



2j c 2mt sinri(7rt) 



log(-t) . (45) 



Then the functions f ^ appearing in the equation f2"T]) 
read 



C(/?) = 



C f (/3) 



i sinh(^) 



2sinh(f + f)sinh(§ - f ) 



(46) 



g 2 ( ? /3 + 27r+f )S 2 (-i(3+% ) 

_J 5f(2^ + f) 

cosh I S 2 (i(3 + 3ir+ *Z)S 2 (-ip + 7r + %) 



These functions satisfy the following important for us 
properties 



509) 



(HP) 



= 5+08), 



a-P) ^' ch-p) 

CO?) = C(2ttz - 0) , C t (/3) = C t (2 7 ri-/3) 
C( y 9)C+(/3 + i7r)=Asinh('|-^ 



(47) 



JV 



2j~l 
1 

N 



5 2 - 2W (2.+ f) 



4sinh(f + ^)sinh(f 
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